Factorization of and Determinant Expressions for the Hypersums 

of Powers of Integers 

Jerome Malenfant 

American Physical Society 
Ridge, NY 
(Dated: January 20, 2013) 

Abstract 

We derive a compact determinant formula for calculating and factorizing the hypersum polyno- 
mials5f)(iV)^En.=i---EnL ,=i("'-L+i)^ expressed in the variable A^(A^ + L + 1). 
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I. INTRODUCTION 



In this article, we will consider the finite sums and hypersums of positive integers raised 
to a positive integer k. Let Sk{N) denote the sum of such integers from 1 to A^: 

N 

S,{N) = J2^'. (1) 

71=1 

The hypersums sI^\n) are then defined recursively as 

N 

S'^\n) = S,{N), Si'^\N) = J2st'\n) for L > 1. (2) 

n=l 

In particular, 

N N 

Si:\N) = J2sS>\n) = J2n = Sr\N), (3) 

n=l n=l 

from which it follows that Sq^\n) = for all L > 1. In the following, we will 

mostly take A; > 0. ^A^) is an [L + k + l)-order polynomial in A^, and is given by the 



formula 



Si'^\N) = J2S{k,q)q\ I ^ + ^ + ^ I (4) 



q=0 



q + L + 1 



where S{k, q) is a Stirling number of the 2nd kind. As an example, we quote from Knuth 
|l| a partial result of Faulhaber for S'g^°^(A^) (which in Knuth's notation is S^^n^): 

S'J^°)(iV) = ^ I 6 A^^^ + 561 A^^^ + ■ ■ ■ + 1021675563656 A^^ + 96598656000 A^ } . 

For A^ = — 1, —2, . . . , —L — 1, the sum over q on the right side of eq.(4) is zero since 
the binomial coefficients are all zero; the sum is also zero for A^ = if /c > 0, since in this 
case S{k, 0) = and all g > binomial coefficients are zero. The polynomial representing 
sI!"\n) for /c > therefore has zeros at these values, and is expressible in the form 

) (AT) = N{N + l){N + 2)---{N + L + l)x Q^''-'^ (N) (5) 

where Q^'^'^^N) is a (A; — l)-order polynomial in A^. 

While the 5*^^^ polynomial can be calculated using eq.(4), this formula does not readily 
lend itself to expression in the factored form of eq.(5). In addition, it is known that these 
polynomials simplify when expressed in the variable A^(A^ + L + 1). The aim of this article 



is to derive an alternative formula for the S*^^^ polynomial in this variable, one in which 
the zeros at = 0, —1, . . . , —L — 1 are factored out and in which the remaining {k — 1)- 
order polynomial is given as a determinant of a relatively simple matrix, one that does not 
explicitly involve Stirling numbers (or any other "complicated" numbers). We do this in 
Theorem 3 in Section IV; in particular, our result for S^^"* is 



^ ^ yiy + 10)(, + 18)(, + 24)(, + 28)(, + 30) x V±1H^. 



where y = N{N + 11). To do this, we need to develop some machinery, so for now we will 
consider only the L = polynomials. These are given by Faulhaber's formula 

SkiN) = T^y :(-!)'- I ^ ^ ^ I 5fc+i-niV", (6) 



k+1 / 

n=l \ 



n 



where 5„ is a Bernoulli number. In the following section, we will derive an alternative ex- 
pression for Sk which seems to be more convenient in extending the formalism to hypersums. 

II. SERIES EXPANSION 

In the sum over n for Sk+i{N) we make the substitution n — m + 1, add and subtract 
terms to make the sum on m go from 1 to N, and then expand (m + 1)^^+^ binomially, with 
the result 

{k + 1)S,{N) = {N + - 1 - E f ^ ^ I S,{N). (7) 

1=0 \ Q J 

Repeating this procedure but now with the replacement n — )■ m — 1, we get: 

^-1 / k + l\ 

{k + l)S,{N)=N'^'-J2\ i-l)'-'S,iN). (8) 



=0 \ Q 



Adding (7) and (8) and dividing by 2 (A; + 1), the odd-g terms in the sum cancel if k is an 
even integer, while for odd k the even-g terms cancel. We therefore get separate recursion 
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relations for the even and odd power sums: 



2p+l 
1 

2p + 2 



{N + + iV2f+l - 1 



{N + + 7V2P+2 _ 1 



p-1 

E 
E 

g=0 



2p+ 1 

2p + 2 
2g+l 



S2,(iV) 



(TV) 



(9a) 
(9b) 



We now define x as 2A^ + 1 and the functions fp{x) and g'p(x) as 



1 



^p(a;j 



2' 

[(^ + l)2p+2 ^ _ i^2p+2] _ 1 + 1) - 1 



We will assume for now that p > 0. Equations (9a) and (9b) then become 
1 



(10a) 
(10b) 



'S'2p+i(a^) = 



2p + l 
1 

2p + 2 



p-i 

9=1 



2p+l 
2g 

2p + 2 
2g + l 



'S'2g(a;) 



(11a) 



(lib) 



(Note that we have absorbed the q — ^ terms in the sums in (9a) and (9b) into /p and g^^ 
respectively.) 

jpix) is an odd function of while g^ix) is even. Further, /p(x) has zeros at x = 0, ±1, 
and gp{x) has double zeros at x = ±1. S2p{x) and S2p^\{x) are therefore odd and even 
functions of respectively, with zeros at these points. These sums can therefore be expanded 
in these functions, which we do as 



S2, = (2p)!j] 



(1) 



■p-q 



r (2g+i)! 



fq-i 



c. 



(2) 



5.^: = (2p+l)!X:p^9,. 



fl2a) 



(12b) 



Equations (11) fix the values C^^^ = c'i^ = 1. We now substitute (12a) into (11a) and (12b) 
into (lib) to get 



(1) 



(2) 



p-1 — p 
^(2g + l)! ^f2a^ 



1^ = (p > 0). 
(2g+l)! ^ 



(13a) 
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These two sets of coefficients satisfy the same recursion relation and have the same initiahzing 
value; they are therefore equal, and we will denote the pth coefficient simply as Cp. These 
coefficients satisfy the additional recursion relations: 

1 



q=0 



1 



(2g + 2)! {2p+l)V 



q=0 



E. 



J2 2'^-'^Cp^q 

q=0 
P 

^2q^p-q _ <2^V(J 

V 



2p 



(2g + l)! {2p)V 



E2qCp-q 



{2q)\ 2i-2p-r 



ECp-q 
(Or 

q=0 
P I 

^ (On 4- 

q=0 

P I 

q=0 
P 



Cn 



a 



p-q 



'Op 



(2g+ 1)! 22P-29 - 2 



p-q 



Cn 



2(2p)! 
1 



{2q)\ 22P-29 - 2 22p - 2 2(2p - 1)! ' 
2p-l 



2p ^P' 



9=0 



and in general, for an integer > 1, 

p 



E 

9=0 



E 

9=0 



a 



p-q 



(2g+l)! 



{N - + (A^ - 3)'*' + ■ ■ ■ + 



2p 



(2g)! 



+ 



(2i-2p _ 
2 



Cn 



{2p-l)\ 



2p-l 



(13b) 
(13c) 
(13d) 
(13e) 
(13f) 

(13g) 
(13h) 
(13i) 
(13j) 







N even 


\ 22P 


1 


N odd {p > 0) 



(13k) 
(131) 





1 


N even 1 






AT odd j 



In fll3ep and fll3fp . -E^ is the kth Euler number. For the derivation of relations (I13bm3ip . 
see Appendix A. 
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Consider now the infinite sum Yl'^=o Cpx"^^, which we assume to be convergent in some 
region. From (13a) we have 



2p 



p=0 



oo p— 1 
p=l q=0 



(2p-2g + 1)! 



X 



2p 



X 



2p 



^ ^ ^ (2)0 -2g+ 1)! 
p=q+i ^ ^ >i I ; 



g=0 71=1 ^ ^ 



g=0 



sinhx 



- 1 



X 



(14) 



from which we get 



p=0 



2p _ ^ 

^ sinh X 



(15) 



As a consequence, the Cp coefficients are related to the Bernoulh numbers as 



2 - 



Bo 



(16) 



{2p)\ 

and correspond to integer sequences A036280 (numerators) and A036281 (denominators) 
in the Encyclopedia of Integer Sequences j^. The Bernoulli numbers are related to the 
Riemann zeta function ({s) for s an even, positive integer 



:|8 



p. 12): 



^^. = (-ir^|^c(2.). 

Combining this with ( !T6|l . 

Cp = (-1)^4^1 -2'"'"K(2P)- 



(17) 



(18) 



But since 



V 2^ / ~ \ 2^ / ^ ~ 2^' ^ 3^' 4^ ^ 5^ 



we have 



1111 

-'■~22p~^32p^42p~^52p 



(19) 



(20) 



where //(x) is the Dirichlet eta function. 



Recursion relation (I13ap corresponds to eq.(3.2) in Van Malderen [6|], and these Cp co- 
efficients are, up to a factor of (— 1)^, the same as Van Malderen's Dp coefficients in that 
article; (see also the article by Chen [7|). Relations ( |13gP and f llSip . for = 3, correspond 
to eqs. (2.2) and (2.3), respectively, in Van Malderen. 

Theorem 1 For k > 0, Sk{x) is given by the expression 



Proof. The sums in eqs. (12) for 5*2^ and S2p+i can be extended down to g = since 
/o = fl'o = 0. The resulting sums are however not valid for the p = sums 5*0 and ^i. They 
can be made so by adding and subtracting terms inside the summation which sum to zero 
for p > as a consequence of recursion relation fll3ap . The "corrected" sums are: 

Expanding the functions inside the brackets in powers of x, we have: 



+ 1)29+1 + _i)2,+i 1 1 ^/2g+l\ 1 




(22) 




and 

{X + 1)29+2 + (x- 1)29+2 1 1 
229+3 

Then, 



a;2"+' - 1). (23) 



(2n+ 1)! ^ (2g- 2n)! ' 

n=0 ^ ' q=n ^ ^ ' 
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The sum over q is evaluated using (113dp : 



{2q-2n)\ ~ (21)1 ~ W^' ^ ^ 

and so we have 

^^f. - 1 _ (2p)! 

' ^ ""^-^ (2n + 1)! A. ^5 (2p + 1 - 25)!' ^^""^ 

Similarly, for 6*2^+1: 



p 



2 ^ (2n + 2)! ^ (2g-2n)! ■ ^ 

The sum over q is the same as in (24) above, and so we get 

(2p+l)! ^ x^"+^-l _ (2p + l)! ^ ^2p+2-2,_i 
^2p+iW- 22P+2 Z.^P-"(2n + 2)! ~ 22P+2 Z.^'? (2p + 2-2g)!- ^ 

QED 

It is interesting to compare the expression in Theorem 1 with Faulhaber's formula. Con- 
verting to Bernoulli numbers and expressed in the variable N, it becomes. 

Sk{N) = ^^^^^ E ( ^ ^ ^ ) (1 - 2'-")i?.+i-„[(2iV + 1)" - 1], (29) 



n=l 



The power-sum polynomials take a simpler form in the variable y = N{N + 1) compared 
to X. (The fact that the polynomials for the odd-power sums simplified when expressed in 
thid variable was known by Faulhaber; see the discussion by Knuth [1].) In this variable, 
the polynomials are: 



_ VW+^ 1 Q _y 

^° - 2 ~r ^^"2' 

^2 - -Q ' ^3 - ^, 

. VWTTyi3y-l) ^ y^2y - 1) 

30 ^ 12 ' ^ ^ 

_ V%+T y{3y^ - 3y + 1) _ y^Sy"" ~ Ay + 2) 

Oh — r:: , O7 



52 ' 24 ' 

_ VW+l yj^y' - lOy' + % - 3) „ _ y\y - l){2y^ - 3y + 3) 

- 90 ' ~ 20 ■ 
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Corollary 1 The power sum of odd integers, 
Sk{N) = 1 + 3^= + 5^= + ■ ■ ■ + AT^ = odd integer. 
is given by the expression, 



k\ 



\k/2\ 
9=0 



(iV + l)*^+l-29 

(A; + 1 -2g)! ' 



Proof : This follows directly from relations fll3kp and (llSip . 
Therefore, expressed in the variable 

NiN + 2) 

y = -^ ' 



(31) 



the Sk are, up to an overall constant and a factor of 2^, given by the same polynomials as 
in (30): 



Sk{y) = 2' 



Sk{y) - Sk 



(32) 



III. MATRIX FORMALISM 



Recursion relation f ll3al) can be written in triangular-matrix form as 



/i 



i i 1 

5! 3! 







C2 



\ 



J 



\ '■ J 






V-/ 



(33) 



Lower-triangular matrices which have I's along their main diagonal are unit lower-triangular 
matrices. Matrices which are constant along all diagonals are Toeplitz matrices. The matrix 
above is thus a unit lower-triangular Toeplitz (LTT) matrix. 

If A is an n X n lower-triangular matrix, {n < oo), then, for k < n, the k — truncation of 
A is the k X k lower-triangular matrix obtained from A by removing all rows and columns 
greater than k. It is straightforward to show that the product of two truncated lower- 
triangular matrices is equal to the truncation of the product and, consequently, the inverse 
of a truncated lower-triangular matrix is the truncation of the inverse matrix. 
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J 



We define the infinite-dimensional lower shift matrix J as 
/ n \ 







1 
1 



(34) 



/ 



and, for future reference, the unit column vectors Ip 



ii 






\''J 



(35) 



Any infinite-dimensional LTT matrix A can be expanded out in powers of J, with J° being 
the identity matrix: 



A 



\ 







q=0 



(36) 



ao 

tti ao 

02 Ol CtO 

V ^ ■■/ 

LTT matrices commute with one another since commutes with J^. The determinant of 
a finite-dimensional lower-triangular matrix is the product of its diagonal elements, so all 
finite unit lower-triangular matrices have determinant 1. 

We will be concerned in the following with determinants that have the general structure: 







Co 

Cl 



(37) 



(lk-2 do Ck-2 

afe-2 • • • oi Ck-i 

To describe them, we introduce some notation: The determinant above is oi a k x k matrix 
consisting of the sum of a "almost LTT" matrix, one whose (k, k) element is zero, (the 'base' 
matrix), and a matrix whose only non-zero elements are along the last column, (the 'tower' 
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matrix). For < q < k — 1, we define the "unit tower matrix" Kg to he a k x k matrix 
whose elements are all zero except for the {k — q, k) element, whose value is 1: 



{k,k-l) 



\ 



\ 







(38) 



(The dimension k of these matrices will usually be clear from the context.) Matrices such as 
A in (36) and others defined below are infinite-dimensional, whereas in taking determinants 
the matrices are finite-dimensional. To keep the notation simple, when taking determinants 
we will use the same symbols for the truncated, finite-dimensional matrices as the infinite- 
dimensional ones, but will indicate with a subscript the dimensionality of the matrices inside 
the brackets. In this notation, the determinant in (37) is 

det |>1 - aoKo + J] Cfc-i-gi^J . (39) 

We can also use this notation when ^4 is a non-Toeplitz, lower-triangular matrix with con- 
stant diagonal element Oq. 

Some useful identities are given in the following lemma: 

Lemma 1 Let A be a lower-triangular matrix which is constant (— Qq) along the main 
diagonal, T be an arbitrary tower matrix, L be a unit LTT matrix, and D be a diagonal 
matrix: D — diag {di,d2, ■ ■ ■}■ Then 



(I) 
(11) 

(III) 
(IV) 



det{A-aoKo + xT}^k) = x det{A- aoKo + T}^k), 



k-l 



k-1 



det <j ^ - aoKo + ^ Ck-i-qKg 

9=0 >> (k) 

det{A - aoKo + LTj^k) 
det{A - aoKo + DT}^k) 



^ Ck-i-q det {A - uoKq + i^J(fe) , 

g=0 

det{L-M-aoiCo + T}(fe) 

4 dei{D-^AD - aoKo + T}^k) 



Proof: The first identity is a statement in this notation of the well-know property that a 
common factor of any column in a matrix can be factored out of its determinant. The 2nd 
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identity corresponds to an expansion by minors of the determinant on the left along the kth 
column. The 3rd and 4th identities follow from the identity (detX)(dety) = det{XY), and 
the easily demonstrated relations ( /c x A; truncation assumed): L~^Kq — Kq, D~^KqD — Kq, 

and TD = d^T. 

In particular, for tower matrices Ti and T2, we have from (II), 



det {A - aoKo + Ti + Ta}^^) = dct {A - a^Ko + Ti}^^^ + det {A - aoKo + T2} 



(fc) 



(40) 



We define the matrices P and Q as 



P = 



Q ^ 



^ 1 



3! 

^ ^ 1 
5! 3! ^ 



h 1 
4! 2! 











\ 



q=0 



(2g + l)!' 



J" 



7 



Then the Cp coefficients are 



(41a) 



(41b) 



Cp^det{P-Ko + Kp} 



(fc) 



(42) 



for any k > p. This is demonstrated by noting that Co = 1 by this formula and that 



E 



a 



p-q 



- E7^;;^det{p-Xo + W 



o(2g + l)! ^(29+1)! 



(fc) 



= det ^ P 



J" 



(43) 



(fe) 



is zero since the last column and the (p + l) column in the determinant are equal. It follows 
from (42) that 



q=0 



(44) 
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Written out, and expanding by minors, the RHS of eq. (42) is the p x p determinant 

1 



3! 

5! 3! 



••• 

1 ••• 



1 



1 



1 



(2p-l)! (2p-3)! (2p-5)! 
Ill 



1_ 

3! 



(45) 



(2p+l)! (2p-l)! (2p-3)! 

which is the form, up to the factor (— 1)*^, that Van Malderen gives for his Dp coefficients 
We now apply this formahsm to the power sums Sk{x) for k > 2: 



^2p[x) - 1^;^ 2^^v~n I 1)1' ^iv^A^) - 22p+2 Z^'-^p-""?: 



2n+l 



(2p+l)! 



2n+2 _ 1 



n=0 



(2n + 



?i=0 



(2n + 2)! 



(46) 



We now define the diagonal matrices as 



ml 



\ 



(m + 2)! 

(m + 4)! 



V 







(47) 



where m takes on nonnegative integer values. These matrices satisfy the identities 
D^'jPh = [(m + 2p)!]±i J^Ii, 



(48) 



as well as the corresponding identities for the truncated matrices, with Ii replaced by an 
appropriate unit tower matrix. 
We have 

(2t)V -J^ ^2n+l 



n=0 



det < P 



Ml det<lP-/^o + y-^-— -p 
22p+i I " 2n+l ! ^ 

(2p)\ X 

22p+l 

(2p)! X 
22p+ 



det <^ p - i^o + ^r^ XI 



n=0 



(P+1) 



(P+1) 



(P+1) 



(49a) 
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And, in a like fashion, 

(2n + 1 V ^ ^2n+2 _ 1 

n=0 ^ ' 



(2p + l)!(x^ - 1) det ^ P - + ^ i^p ^ . (49b) 



Note that 



(p+i) 



(7 - J)DiPD^Hi = {I - JfD2PD^% = Ii (50) 

I.e., the first column of (7 — J)DiPD^^ and of (7 — J)^7)2-P-DJ^ consists of O's except for 
the (1, 1) element, which is 1. As a consequence, 

det \p-Ko + D^' -1- Kp] = det {(7 - J)D,PD:[' -Ko + Kp} = 0, 
I ^ 'J ) (p+i) 

(51a) 

det |p - TTo + D^' -^f^j^ ^p} = det {(/ - JfD^PD^' -Ko + = 0. 

(51b) 

Using these identities, we can thus subtract off the tower matrices D'^^{I — J)~^Kp and 
75^^(7 — J)~'^Kp from the tower matrices in (49a) and (49b), respectively, without changing 
the values of the determinants, to get. 



(52a) 
(52b) 



The matrix J, acting on Kp inside the brackets, lowers it to Kp^i. The (l,p+ 1) element of 
the matrix in each determinant is therefore zero, and we can expand each by minors along 
the first row, reducing them to p x p determinants. Changing to the variable y and using 
identities (III) and (IV) in Lemma 1, we get: 



^2p(^) = 2^v-^^i) det {[ / - (4y + 1) J ](7 - J)D,PD^' -Ko + Kp_,}^^ ; (53a) 
'^2p+i(^) = ^2p-2(2p + 2) ^ - + 1)-^ - JfD.PD-,' -Ko + Kp_,}^^ . (53b) 

Theorem 3 in the next section generalizes these expressions to the L > cases. 
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IV. HYPERSUMS 



We generalize relation (44) and define the coefficients Cp^\ L > 0, by 

oo 

p-(^+i) = ^Cf)j^ (54) 

p=0 

We then have 

Cr = ^P, ^7(^^ = 1, C7W = f^Cf-i)C,_„ L>1, (55) 
and 

) = det {P^+i - i^o + . (56) 

If we multiply this determinant on the left by det(P'^^^^^)(= 1), where R is an arbitrary 
integer, and take the matrix p^-^-i inside, (56) becomes, with the replacement p — k—l—q, 

Cf\-, = det {P^ -Ko + P^-^-V^i^._i}(,) . (57) 

This second form is more useful in the proof of Theorem 3 below, where we will be adding 
terms containing Cp^^^s with different L and p values ; we can sum these up into a single 
determinant if the base matrices are the same. 

We state below some recursion relations that these coefficients satisfy: 

y \ ^P-^ = y = Ci^)- (58a) 

^ (2g + 2)! ^r^(2g + l)! ^ ' ^ ' 

^(2,+ !)! - (2,)! - L + l • 

(The ffist equality in each line also holds for L = —1.) Using the third equation in (55), the 
two equalities in fl58ap follow from f ll3cp and (llSal) . respectively, and the ffist equality in 
( I58bll follows from (113bl) . To prove the second equality, it follows from (54) that the Cp^^^s 
are the coefficients in the expansion of (x/ sinhx)^"*"^: 

L+l °° 

= 5^Cf)x^^ (59) 

p=0 



X 



sinh X 

Differentiating both sides 



iL + l)x'^ _ {L + l)x^^'coshx _ ^ ML).p-i /gO^ 



sinh X sinh x 

p=0 
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Rearranging terms and expanding out the cosh x factor, we get 

L+1-2P _ x^+^cosh(x) _ ^ x^'^^^ 

h ^ + 1 " sinh^+^x (2?)'' 

The proof then follows by equating powers of x. 
We define the variables xi^ and yi,: 

xl = 27V + L + 1; = 7V(7V + L + 1); ^ = 47/^ + (L + 1)^. (62) 

(To keep the notation simple, we will in general suppress the subscript L on these variables: 
— )■ a;, Dl ^ y, with the L understood.) The generalization of Theorem 1 to hypersums 

is: 

Theorem 2 For L > and k > 0, the hypersum polynomials are given by: 

, p+M f 2o+l 1 ^(2.U+2-2,s) 

c.(2M). X ^ {M>y- ^ J^(2M) 1 ^ S+.u^i-,j-,s ■'-(■'• + ^'^ - -3)! 

-^2^ W 22P+1+2M (2^ + 1)! (2g+l)!^ (2s -1)! (x-2s + l)!! 

C.(2M+1). X ^ (2p)! l^(2M+l) 3^^^ 1 ^p+M+l-q-s X{X + 2s - 2)\ 

^2p W 22P+2+2M |S+M-g ^2g + 2)! (2g+l)!^ (2s)! (x-2s)!! 

(2M-1). X ^ (2p+l)! f (2M-1) _ J_ ^ Cg^^^Vi ix + 2s-2y.\ \ 

^2p+i [X) 22f+i+2M ^2g + l)! (2g)!^ (2s -1)! {x-2s)U J' 

(,M). X ^ (2p+l)! / „(,M) x^^-^^ - 1 _J_^ C%l-:!l^_, ix + 2s-m \ 
^2p+i [X) 22P+2+2M |S+M-g + 2)! (2g)! (2s)! (x - 2s - 1)!! J ' 



where we use the convention that Cp — for p < 0. See Appendix B for the proof. 

The sums over q in the second terms in these expressions can be evaluated using some 

of the relations in (58), but the unsummed expressions above are more convenient in the 
proof of Theorem 3 below. Before getting to that theorem, we first need to prove a technical 
lemma: 
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Lemma 2 The matrix P satisfies the identities, 

(I) D^P^'^+^D^^J'' Ii = 

(II) D2P''+'^2"'^' Ii = ,,w, " n— ^TTTT 





-J){I- 


-9J)---[7 




T 

(2A; + 1)V] ' 












{I 




- 16 J) • • • 


[I 


- {2k + 2)2 J] 












{I 


-AJ){I 


- 16 J) • • • 


[I 


- (2A; + 2)2 J] 










. Ii 



(III) D^QP^'^+^D^^J'- Ii = 7-;^ — „ , „ II, 

(IV) i,.QP-«D,V' I. - _ ^jj^ _ 1^ _ ^ 3j,^j 

Proof. (We remind the reader that the matrices are defined by equation (47).) Each 

equation above is a statement that the first column of the matrix on the left is equal to the 
first column of the matrix on the right. The matrices themselves are not equal since the 
ones on the right are LTT and those on the left are not. 
Let the matrices and be defined as: 
J'' 

= {i-j){i-9j)---[i-{2k + iyjy ^^^""^ 
J'' 

^ (7 - 4 J) (7 - 16 J) ■■■[I-{2k + 2)2 J] ■ ^^^^^ 
These matrices can be expanded, respectively, in powers of odd and even integers as 

k (fc) oo k 

q=0 ^ ^ ' n=0 g=0 



= E ,_(2;+2)2j - E E/^f + 2)^"- (64b) 



q=0 ^ ^ ' n=0 ?=0 

The proofs of the first equalities in (64) follow by induction on A;, and then the inverse 
matrices are expanded in positive powers of J. It is easy to show that these sets of coefficients 
each satisfy [k + 1) independent equations, 

E:=o4'\2g+i)2^ 



1, forp = A;; ^ ^ 

(65) 

0, for < p < k\ 



(project the far left and the far right sides of equations (64) onto Ii: in powers of J, 

and Yfc have the form + a\J^^^ + a2J^^'^ + " ■ ■ ) so the resulting column vector's first k 
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elements are 0, with the {k + 1) elements equal to 1). The expansions for Xk and Yk are 
then: 

oo k 

Xk = ^''^(^g + l)2"+2'=; (66a) 

71=0 g=0 
oo k 

n=0 q=0 

Relations ( |65ll are sufficient to determine the [k + 1) values of these coefficients. The coef- 
ficients are in fact given by exphcit formulas, but for the purpose of this proof, only these 
relations are required. 

We consider relation (I) ffist. We will use the expansion for P in powers of J given in 
(41) and the m = 1 identity in (48). For k = 0, the left-hand side of relation I reduces to 

D,PD^' Ii = Y, J'h = Ii, (67) 

9=0 

which equals XqIi as required. For A; > 0, an evaluation of the left-hand side of the relation 
results in 

D^p^k+ij^-ijk _ ^^j^ (gg) 
where 

_ >^ 7^=+" >A / 2n + 2k+l 

Zk is an LTT matrix. Relation I therefore holds if and only if this matrix is equal to the 
matrix X^ for all k. We will assume then that, for some k, is given by the sum 

oo k 

Zk = J2 J""^^ "i'^^2g + 1)'"+'^ (70) 

n=0 q=0 

Then 

_ {2k + 1)! ^ ( 2n + 2k + 3 



X 

p'^r'o \ 2pi + 2A; + 1 
18 



q=0 



Using the binomial identities 

n 



fc=0 



2k + 1 I 2 



E 

fc=0 



2r2 + 1 \ _ + 1)'"^' -{x- 1)'"^' 



2A; 



and relations dHS]) for the a's, this reduces to 

oo fc+1 

n=0 (j=0 

(2A: + 1)! r 2g + l (fe) 2q + 1 (,) \ 

4(2fc + 3)! \ 2g- 1 "^"^ ^ 2g + 3 ^+^J ^ ^ 

(a^^) = for g < or > k). 

It is a straightforward calculation to show that the coefficients in the sum over q satisfy 
relations f l65|) for, and thus must be equal to, aq''^^\ As a consequence, Zk+i = X^+i and, 
by induction on k, for all k' > k. 

The proof of relation (II) follows by a similar argument, using the matrices D2 and Y^. 

To prove (III), we use relation (I) to write the left-hand side of (III) as 

D.QP^'+'D-'J'' Ii = D,QD-'Xk Ii. (73) 
Evaluating this using the expansion for Q in (41b), we have: 

OD jp 00 k 

p=0 ^ m=0 q=0 

00 k (fc) 

n=0 g=0 ^ ~'~ 

= f; J-'.^(2, + 2)-«{|-ti af + a« |l. (74) 

n=0 q=0 I. y "T y "T J 

Again, it can be shown that the coefficients in this equation satisfy (1651) for /3q''\ and we 
therefore have 

D^QP^^+^D^\J^ Ii = nil, (75) 
which is relation (III). 
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The right-hand side of (IV) is 'almost' Xfc+ili, but it lacks one power of J. Nevertheless, 
we can write the left-hand side, using (II), as 

D^QP^^+^D^^J^ Ii = D^QD^^Yk h, (76) 
which evaluates as 



JP ^ 

p=0 ^ m=0 (j=0 

00 k r,(k) 

= E E ofo\ 9^2 { (2? + 3)^"+^^^^^ + (2g + l?-^''^' - 2} Ii 

n=0 g=0^^^^+^^ 
00 fc+1 r Mk) oik) ^ 

where 

The coefficients in the sum over q satisfy (165|) for a^'*'^ and therefore 

D2QP^^+^D:,^J^ Ii = E E «f ^'H2g + i)2n+2/=+2j^^ (78) 

n=0 9=0 

which is relation (IV). 
QED 

We will use this Lemma in a slightly different form in proving Theorem 3 below. For 
{k X /i;)-truncated matrices J, P, etc, the identities in the Lemma also hold if the unit column 
vector Ii is replaced by the x unit tower matrix Kk_i; both Ii and K^-i project out the 
first column of any matrix that they are multiplied by from the left. We will refer to both 
forms (the "column- vector" and the "tower-matrix" forms) simply as "Lemma 2" . 

Our main result is: 

Theorem 3 Expressed in the variable y = N{N + L + 1), the hypersum polynomials S\^\y) 
for L > 0, k > are given by 

Li/2J 

{L + k + l)\ 2'=-i(l + L mod 2~ 



S^'^(v) fc! (V4,+ (L + 1)2)^ Ti\v + a(L + l-a) 



q=0 

LL/2J 

det <j [ / - (% + (L + 1)2) J ] J] [ / _ (L + 1 - 2qfj ] D^P'^+^Dr^^ -K^ + K^, 

q=0 
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where 

A = L mod 2 + (A; + 1) mod 2 = 0, 1, 2; 
m — L + 3 — k mod 2; 



Proof. We consider first the 5'2p^^' polynomials. Then A = 1, m = 2M + 3, and n = p. 
The sum over q from to p + M in the expression for this polynomial in Thm II requires 
that the size of the matrix used to represent the coefficients be at least {p + M + 1). Using 
(57) with R — 2M + 1, we can express the C^^_g and the C^_^^i^gl.g coefficients as 
determinants with the same base matrix: 

Cm-. - det -Ko + J^K^^m} ; (79a) 

Cm^i-;!. - det {P^--^^ -K, + P''-'J'^'-'Kp^m},^^^^,, . (79b) 

Then we have from Thm 2 and identity (II) in Lemma 1 that 

, p+M ( 2o+l 1 M ^(2M+2-2s) 

(2M). X ^ (2p)! I ^(2M) a:'^^ 1 (^p+M+i-q-s x{x + 2s-3)U 

^2p [X) - 22P+1+2M jS+M-g +1)! (2g + 1)! (2s - 1)! (x - 2s + 1)!! 

- ^WTT)!^ (2.-1)! (X-2.+ 1)!! 

This expression will be reduced and simplified by making a series of replacements for terms 
and factors in the tower-matrix part which are valid inside the brackets of the determinant, 
(but not generally outside of them) . The first replacement is for the sum over q in the 4th 
term: J2qJ'^/i'^Q + 1)! — > P. Although P is given by an infinite sum, this replacement is 
allowed since P inside the brackets is truncated. Other factorials in the denominators can 
be dealt with by making the additional replacements: 



p+M, 



Kp+M D^^J^ ^Kp^M- 



(2.-1)! 
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Eq.(80) then becomes 



-£..P-.rV.-f±|^)A,,„^ . (81, 



(p+M+1) 



The ratio of double factorials is equal to x for s = 1 and is 

4- 9« _ ^ ^^^2 _ ^^^^2 _ 9) . . . (^2 _ - 3)2) (82) 



{x-2s + l)\\ 
for s > 1. We can also make the replacement 

TS-l 

7-) p2s-l 7-)-l TS-1 , ^1 

' ^ (/- J)(/-9J)---[/-(2s-l)2j]- 

This is allowed by the "tower- matrix" form of relation (I) in Lemma 2. Upon making these 
changes, we get 



s=2 



ll/_(2g+l)2J 



1. 



P+M+1 



The terms inside the brackets can be summed as 



X X X ^ 'yix''-{2q-lf ^^x2-(2g + l)2 



7-J 7-jZ-.ll/-(2o + 1)2J 11 / _ (2o + 1)2J 7 - 



,=2 ,=1 ^ - (2? + 1) V / - (2g + 1)2 J 7 - 

(84) 



by successively combining terms with the same power of J: 

X X x{x'^ — 1) J 

J-a;2j ~ 7-J ~ I - J I-x^J' 

x{x^-l) J x{x^-l)J x{x^ -l){x^ -9) J2 



I- J I-x^J (/-J)(/-9J) (/-J)(/-9J) I-x^J' 
etc., 

Inserting (84) into eq. (83), factoring out the a:;n(^^ ~ (2? + 1)^) (as per identity (I) in 
Lemma 1), and using J^Kp^M — Kp (inside the brackets), we have 

{ 1 1 1 



(p+M+1) 
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Now consider the value of the determinant in this equation ai x — 2M + 1. By relation 
(I) in Lemma 2, this is equal to 

This determinant however is zero, since it contains two columns, (the (M+l)-th and the last 
column), which are proportional to one another. We can therefore make the replacement 
1^1 1 _ x'^- (2M + 1)2 J 



7-^2 J I-x'^J 7-(2M+l)2J I -(2M + iyj I -x^J 
in the determinant in (85) without changing its value. Factoring out x'^ — (2M + 1)^, this 
equation then becomes: 

q=0 

This expression accomplishes, for these values of L and k, the goal set out in the Intro- 
duction, to derive a determinant formula for the hypersum polynomials which factors out 
the zeros at A?" = 0, . . . , —L — 1. The size of the determinant can however be reduced; the 
first (M + 1) elements in the last column of the tower matrix in (87) are zero and we can 

make the replacement in that equation of 

0, 



'(M+l) 



(88) 



Successively expanding the determinant by minors along the top row, (87) is thus reduced 
to the p X p determinant: 

M 

'2p \-^J — 22p+l+2M 



4rw = ^Mif-ni-'-p'+i)'] (83) 



5=0 

X det 



Now applying identites (III) and (IV) in Lemma 1 and changing the variable to y using 
x'^ — Ay + {L + ly , eq. (89) takes the form as stated in the Theorem for L = 2M, k = 2p. 
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Likewise, the expressions in Theorem 2 for the other hypersums can be written in deter- 
minant form as: 



M 



(p+M+l) 



c(2A^)^^N _ (2P+ Het / p2M+i _ 7^ , r)-i 

'^2p+l \X) - 22P+2+2M Clet < Ao + 1^2 |^(J _ J)(J _ ^2 J) 



-A/ 



(p+M+1) 



Using relations (II), (III), and (IV), respectively, in Lemma 2, they can be brought into the 
forms stated in the Theorem by a similar process. For S'^p_ll and •Sgp+j"^'*, the determinants 
corresponding to that in (86) for 5*2^^^'' are not equal to zero, and so the final forms for these 
polynomials contain the tower matrix Kp rather than Kp_i. {q] 
QED. 

We have already given the expression for S^^^^ in the Introduction; some other examples 



are: 




sf\y) = 


6! 
10! 


Sf\y) = 


7! 
13! 


Sl?{x) = 


11! 
20! 


sit\y) = 


14! 
29! 



— . , 7y^ + 14y=^- 238|/ + 295 

, 14t/5-4011?/3 + 25868?/2 + 145896?/- 1199616 
y + 8)iy + U)iy + 18)(y + 20) x ^ ^ ^ 

4u + 225 

^ yiy + 14) (y + 26) (y + 36) (y + 44) (y + 50) (y + 54) (y + 56) 

X {y^- 1750y^ + 29960?/^ + 376167?/^ - 11436860?/ + 62455917 ) 

Now let A^^'' be the determinant in the expression for 5f ^: 

{LL/2J 
[ / - (42/ + (L + I f) J ] J] [ / - (L + 1 - 2g)V ] D^P'^^'D;^' - Ko + K^^-i 
> in) 

where m and n are as defined in Theorem 3. 
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Proposition 1 a1^^ has the series expansion 



s=l 



where a, is the s x s determinant 



LL/2J 



as = det {[I-(L + 1) J] [7 - (L + 1 - 2g) V] L'L+fe+3-2.P''+'i?Zife+3-2. - + K^-i 

q=0 



Proof : We re-express A\^^ as 



LL/2J 



g=0 



[4y + {L+iy]J 



Kn-l 



(n) 

(90) 



We have 



1 



oo g 



^[4y +{L + 1)^] = E E r I + If'-'' J', (91) 



/-[% + (L + l)2]J 



g=0 g=0 s=0 \ -S 

although only the q < n — 1 terms in the g-sum will contribute to the determinant. We 
interchange the order of the sums and take the sum over s outside of the determinant: 



n-l 



LL/2J 



A^Hl/) = Y^i^yy det <J n [ ^ - + 1 - 2?)''^ ]DmP'^'''D-' - Ko 

q=0 



s=Q 



n-l 



q=s 

Using the negative binomial series, 

r=0 \ * 



(92) 



(n) 



1-Z 

the sum over q becomes 



(93) 



n-l 



E 



q=s 



s 



i-{L + iyj 



5+1 



(94) 
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We now make the change of summation index s — )■ n — s and get 

n 

Af^(z/) = E(42/)""^«- (95) 

s=l 

where 

as = det l[I-{L + 1) \{[I-{L + l- 2qfj]D^P'^+^D-} - Ko + K^^i 

[ 5=0 

As was done in the proof of Theorem 3, the size of the determinant can be reduced from 
n X n to s X s by the replacements — )■ -Dm+2n-2s- We recall from Theorem 3 that 
m = L + 3 — k mod 2 and n = [{k + 1)/2J ; then m + 2?7, — 2s is equal to L + k + 3 — 2s and 
the result follows. 
QED 

V. CONCLUSION 

We have derived a formula for the hypersum S"^^^ polynomial in a factored form which does 
not explicitly involve either Bernoulli or Stirling numbers. The hard part of the calculation 
has been reduced to the multiplication of matrices constructed entirely from 'simple' numbers 
and to the calculation of a determinant of size \_{k + l)/2j. 
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Appendix A: Recursion relations 

As noted above, recursion relation ( I13ap follows directly from inserting either (12a) into 
(11a) or (12b) into (lib). The demonstration of most of the remaining relations is more 
transparent (or at least more compact) if these are written in matrix form. Matrices P, Q, 
and J and the column vectors Ip were previously defined in section III; we also define the 
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diagonal matrix, 

2^ 

2' 

2« 



2 



(Al) 



and the column vectors, 



C = ^^CgIg+l, C = ^^CgIg+l, Cg 

g=0 g=0 



2g 



(A2) 



{2q)V 

The matrix elements of are proportional to the Bernoulli numbers and those of 



are proportional to the Euler numbers: 

Q-' -W-H^.J'- 

In this matrix notation, f ll3ap is 

c = p-^ii. 



(A3) 
(A4) 



(A5) 



Now consider the product PQ: 



JP+q 



,,.o(2P+l)'(2^)! ^(2n-2g + l)!(2g)! 

- y-_^y- ( 2^ + 1 
~ ^ (2n + 1)! ^ I 

n=0 ^ ' q=0 \ -^Q 

The evaluation of the sum over q is straightforward and equals 2'^"', so we have 



n=0 



In similar fashion, the products and can be evaluated: 

2'in+l jn 



n=0 



{2n + 2) 



n=0 



2'2n+ljn+l 

(2n + 2)! 



I + JP\ 



(A6) 



(AT) 



(A8a) 
(A8b) 
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On the other hand, 

2J''2-^ = 2^T, (A9) 
and therefore 

°° tj) jn f) — l 
n=0 ^ ' 

Also, 

^ (2n + 2)! ^ (2p)! 2 \^ (2p)! /' 

or 

2Q2-^ = I + 2- JP\ (=2-g2-/from (A8b)). (A12) 
Now relation ( ]13b[) 

corresponds to the matrix equation 

^P^^^C = gii (A14) 

and this is derived by multiplying both sides of the equation C = P^^Ii on the left by 
PQ{= QP) and using relation (AlO). The third relation, f ll3cp . 

^(2g + 2)! {2p+l)\ ^ ^ 

is equivalent, under the substitutions p = p + 1 and q = q + 1, (and then dropping the 'bar' 
notation), to 

E 2^^Cp-g ^ ^ 2 
In this form, f ll3cl) is (A12) above multiplied on the right by C: 

2Q2-^C = C + 2 ■ JPli. (A16) 
The next three recursion relations 
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P o2p— 2(7/^ 771 

^ (2g+l)! - (2p)!' 



q=0 
P 



g=0 



all follow from relation (AlO). In matrix notation, fll3dp is 
{2-^Q2)C = 2"^C, or Q2Q, = Q,. 

( ]13dl) is proven by noting that the left-hand side of the equation on the right is 

Q2Q = {p-^2P2-^)2{p-%) = p-^21i = p-% = C, 

(since 2Ii = Ii). Now inverting and rearranging equation (AlO), we have 

2P-^2-^P = Q-^. 

The left-hand side of this equation is 

2P'' 2-'p = yry 2^p-^ \ ^"-\^ , 

^ ^ 2g+l!' 

p=0 g=0 \ ^ ' ) 

while the right-hand side is given by (A4). Equating powers of J, we have 



^P- 

„ (2g+l)! 



<?=o 



which is relation fll3ep . Rearranging the factors in (A20), it becomes 

2-^Q-^p-^2 = p-\ 

From (A3) and (A4), this is relation ( ]13fll . apart from an overall factor of 2" 

p 

Z LjqLyp-q — Op. 

q=0 



2p. 



To prove relation (|13gP 



P ^ 

(2q)\ 21-2P-1' 
consider the matrix 
M = 2- P2-^P-^Q-Q 
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The right-hand side can be rewritten as below: 

M = 2- {2-^PQ)p-^Q - Q, (using P2-^ = 2^^QP) 

= 2-\2Q2-^ + I)-Q, (using 2 ■ = 2Q2-^ + /) 
= iQ + I)2-'-Q. 

Therefore, 

A^li = {Q + 1)2-% - Qh = Qh + li - Qh = li; 



or 



Ii = (2- P2-'P-'Q-Q)h. 

Muhiplying both sides of this equation by P^^, it becomes 
C = i2-2-^-I)QC. 



Solving for QC, we get 
1 



QC 



which is ( |13gp . We now write the inverse matrix in this equat 



2-2-^-1 ~ 2-2-r 
Equation (A31) then becomes 

1 C = --(C + QC). 

2-2-1 2^ ^ V ; 



Recursion relations ( ]13hl) 



^ 1 
^ (2q+1)\ 



a 



p-q 



(2q 
and ffT3il) . 



22p-2g _ 2 



2{2py: 



p-q 



Cr, 



q=0 



{2q)\ 22p-2g - 2 22p - 2 2{2p -1)1' 
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are obtained by multiplying this equation by P and by Q. Multiplying by P we get 



Pj^^C = -^(^C + PQC) = -^(/ + Q)li 



which is (I13hp . If we multiply (A33) instead by Q, we have, 



Q{2 - 2 ■ I)-'C = -^{QC + Q^C). 

From (A8b), Q^C = {I + JP^)C = C + JPIi, so 

Q{2-2- I)-'C = -l-iQC + C + JPIi) = {2-2- ly'C - ^ JPIi, 



(A36) 



(A37) 



(A38) 



which is fll3ip . 

By eq. (54), P^^Ii = C*^^\ where C*^^^ is the column vector whose elements are the Cp^-* 
coefficients. Therefore, multiplying Eq.(A29) by P~^, we get 



(A39) 



or 



gc(^) = 2- {p-^2-^Q)p-%-c^^^ 

= 2 ■ {2-^p-^)h - C^^) (using P~^2-^Q = 2-^p-^ from (A20)) 
= (2-^-i-/)C« 

which gives us 
p M^) 

q=0 ^ 



Comparing this relation to the second equality in (I58bl) for L = 0, we get 



(,) _ 2p-l 



which is (13j ). 

To prove relations ( 113k[) and ( 11311) . consider the polynomial 



E 

q=0 



a 



P^ ^2,+ l_ 



(2g + l)! 



Using relation f lT6|) . this is expressible in terms of Bernoulli polynomials, 

p 



^p-1 ,,29+1 



g=0 



(2g + l)! 



(2p + l)! 



B2p+l{u) — 2^^i?2p+l 



2/J 



(A40) 



(A41) 



(A42) 



(A43) 



(A44) 
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As a consequence of the identity Bk{u + 1) = Bk{u) + ku ^, the RHS becomes 

^(2g + l)! (2p)! ^ ^ A.(2g+l)!^ ^ 

^ {(u - + - 3)2p + . . . + (m - 2Ar + 



(2p)! 



Differentiating this equation with respect to u, 
^(2g)! (2p-l)!^ ^ ^(2?)' 



(2p 



Setting u = 2N and u = 2A^ + 1 in these two sets of equations prove, respectively, relation 
( ]13k[) and, using p3gP , relation ( 11311) . 



Appendix B: Proof of Theorem 2 

The proof is by induction on L and uses the expressions for the integer power sums from 
Theorem 1, 

p 



s (N) - (2^)' i^N+^r'^'-^ 



22p+i^ (2g+l)! 

. (Mf (fap>0), (Bla) 

22p+l P 9 (2g + 1)! ^ ^' ^ ^ 

(2p + l)! ^ (2Ar + 1)^^+^ - 1 

and the odd- integer power sums from Corollary 1, 

^2p(A^) - -^2-.^P~^ (2g+l)! ' (^^"^ 
52p+i(A^) = ^ Z-^^P-g (2g + 2)! ' ^ ^ 
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where in the last two equations N is an odd integer. The expressions for Sjf'^ in Theorem 
2 reduce to the integer power sums above when L — 0. It is necessary then to show that 
these expressions satisfy 



N 



(B2) 



n=l 



Since Sl^\n) = for n = 0, —1, —2, . . . , —L — 1, we can make the replacement 



N 



N 



(2M,2M+1), 



[n) ^ 

n=l n=—M 



E si 



(2M,2M+1) 



n] 



without affecting the values of the sums. Taking the sum from n 
first term in the expression for S^^^ in Theorem 2, we have 



(B3) 

-M to n = N oi the 



N p+M 

E Eq 

n=—M 0=0 



(2M) 



X 



2q+l 



p+M-g (2g + l)! 



p+M ^(2M) 



x=2n+l+2M q=Q 



- E(^^ W2iV+l + 2M) 



q=0 s=0 
^ p+M ^2s+2 P+M 



j^2s+2 



(2s + 2)! 



= E 



E^(2M) ^ 
^p+M-q'^q-s 



2 ^(2. + 2)! 

p+M 

- Vr^ 
9 2^ 



y(2M+1) 



j^2s+2 



2 ^ (2s + 2)! 



X=2N+2+2M 
X=2N+2+2M 

(B4a) 



X=2iV+2+2M 



In the 2nd line we've used (Bid) for the sum of odd integers raised to an odd power. By 
similar calculations, we get for the other sums, 



N p+M 



EST r^2M+i) 

%=-M q=0 



X 



2q+2 



'p+M-q 



(2g + 2)! 



p+M+l 



x=2n+2+2M 



N p+M 2a+l 

ESr^ ^(2M-l) X 
l^^p+M-q (2g+l)! 

29+2 _ ^ 



n=-M+l q=0 
N p+M 

E EC, 

n=—M q=0 



x=2n+2M 



(2M) 



X 



1^ P+M-q {2q + 2)\ 



x=2n+l+2M 



i E o, 

p+M 

lE^. 

p+M+l 

\ E 



^2s+l _ X 



{2M+2) 

P+M+l-s (2S + 1) 



(B4b) 



X=2W+3+2M 



(2M) 



^2.+2 



P+M-s (2s + 2)! 



(B4c) 



X=2Af+l+2M 



v2s+l r)2s V 
(2M+1) ^ — Z A 



s=0 



P+M+1-. (2s + 1)! 



(B4d) 

X=2JV+2+2M 



We now consider the 2nd terms in the expressions for S\^\ We can use the identities 



ix + 2s- 1)!! 



{x-2s-l)\\ 
(x + 2s-2)!! 



2^^ (2s)! 



{x-2s)\\ 



x=2n+l+2M 



x=2n+2+2M 



2s 



= 2 



2s-l 



(2s-l)! 



n + M + s 
2s - 1 
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N 
j=s 

and 



J 
s 




n + M + q + a 
2q + a 

to derive the relations: 
^ ,t(,t + 2s-3)!! 



2(7 + 1 + (/ 
2n + 2M + l + a 



n + M + q + l + a 
2q + l + a 



N 

E 



n=— M 
AT 

E 

n=-M+l 



(a;-2s + l)!! 

x(x + 2s-2)!! 
{x-2sy.\ 



1 X(X + 2s-2)!! 



a;=2n+l+2M 



2(2s) 



{X-2s)\\ 



X=2N+2+2M 



X{X + 2s-l)\\ 



x=2n+2+2M 

x + 2s-2y.\ 



(x - 2s)!! 



2(2s + l) (X- 2s -1)1! 

1 (X + 2S-1)!! 



X=2N+3+2M 



x=2n+2M 



E 



a; + 2s-l)!! 



n=-M 



(a;-2s- 1)!! 



2(2s) (X-2S-1)!! 
1 (X + 2s)!! 



n + M + g + a 
2q + a 



(B5a) 
(B5b) 
(B5c) 



X=2N+1+2M 



x=2n+\+2M 



(B5d) 



X=2N+2+2M 



2(2s + l) (X-2S-2)!! 

With the sets of relations (B4) and (B6) it is straightforward to show that the expressions 
in Theorem 2 for the hypersums satisfy condition (B2) above. For example, with x = 
2n + 1 + 2M and X = 2A^ + 2 + 2M, 

N N 



n=l 



n=-M 

(2p+l)! 

22p+2+2M 

(2p+l)! 



N p+M 



„2g+2 



n=—M 9=0 

p+M+1 

22p+3+2M 

(2^+1)!^+^-^^ 



^P+M-g (2^ + 2)1 (2g)! 



-y 



M ^(2M-2s) 



p+M+i-g-s (a; + 2s - 1)!! 



s=l 



(2s)! (x- 2s -1)1! 



,(2M+ 

p+M+i-g (2g+l)! (2g)!-fr^ (2s + 1)! (X-2S-2)!! 



C. 



(2M+1) - 2^^X 



M ^(2M-2s) 



p+M+l-g- 



(X + 2s)!! 



E 

9=0 

p+Af+1 

22p+3+2M 

9=0 



22p+3+2M 

(2p+ 1)! 



(2M+1) 



J^2</+1 



^ V r£ 



M ^(2Af-2s) 



+M+l-q-s 



{X + 2s)\\ 



p+M+i-9 ^2g+ 1)! (2g)!^ (2s + l)! (X-2S-2)!! 



(2M+1) 



M+1 ^(2iU+2-2s) ... 



p+M+i-g + 1)! (2g)! ^ (2s - 1)1 (X - 2s)!! 



_ q{2M+l) 

— ^2p+l 



(B6) 



where we've used the relation 



^(2^ + 1)' 



p ML+l) 



(B7) 
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from relations f l58bp to move the term in the first summation to the second summation 

where it becomes the s = term, and then relabeled the index, s — t- s — 1. 

QED 
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